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This remark is part of an ongoing project to simplify the structure of the multi-loop anoma- 
ly I lous dimensions for parton distributions and fragmentation functions. It answers the call for 
a "structural explanation" of a "very suggestive" relation found by Moch, Vermaseren and 
■ Vogt in the context of the a; — > 1 behaviour of three- loop DIS anomalous dimensions. It also 

K*" ! highlights further structure that remains to be fully explained. 
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Abstract 



1 Introduction 



This letter stems from a project to better understand the structure of multi-loop anomalous 
^ dimensions both for parton distributions and fragmentation functions [1]. These distributions, 

Q_i' which we shall generically denote as D [Dn{Q'^) in moment space, D{x, Q^) in x-space) satisfy 

^ , a renormalisation group equation 

(D ! 2 

^^^^d,D = ^{N,as{Q'))D^{Q'), (1) 

•rH , 

^ I where dt is a compact notation for the derivative with respect to t = InQ^ and 7(A^, as) are 

' elements of an anomalous dimension matrix. The latter have been calculated in terms of an 

expansion in the coupling Os up to three (two) loops in the space-like (time-like) cases [2-5]. 
They become increasingly cumbersome beyond leading order. 

Conventionally one defines parton splitting functions, P{x), as the inverse Mellin transform 
of the corresponding anomalous dimensions, giving evolution equation in x space in terms of 
a direct convolution 

dtD{x, Q2) = [' ^P{z, as{Q')) D (-, Q2) , (2) 

where D{x,Q'^) has the physical support x < 1. We have reason to suspect that there might 
exist a reformulation of the evolution equations © in which, by generalising the structure on 
the right-hand-side, one is able to simplify the splitting functions. This is equivalent to stating 
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that the higher-loop structure of the anomalous dimensions in can in part be understood 
as inherited from non- linear combinations of lower loops. 

In this picture, the new splitting functions would not only be more compact, but they 
would also exhibit some important physical properties: beyond first loop they should vanish 
at large x, and they should be identical for space-like and time-like evolution, thus restoring 
Gribov-Lipatov reciprocity [6] (broken beyond first loop in the standard formulation, see [4]). 

A possible reformulation of is^ 



This is precisely the relation that was noted by Curci, Furmanski and Petronzio in [4] for 
non-singlet quark evolution. In the singlet case (both for quarks and gluons) there are also 
interesting patterns, see appendix. Furthermore, if one writes V{z,as) as a series in the 
physical coupling, Oph = Oms + ^c^ms + ' ' ' (with K = — ^)Ca ~ ^iT-f): the two- loop part 
of 'P{x, as) vanishes as (1 — x) for x — > 1. This corresponds to the wisdom of Low, Barnett and 
Kroll [9] according to which the classical nature of soft radiation reveals itself at the level of 
the 1/(1 — x) and constant terms. This classical nature of soft radiation allows one to absorb 
all soft singularities into the first loop and to look upon higher-loop splitting functions as due 
to true multi-parton (quantum) fluctuations. 

Finally, one notes also the z-dependence of the argument of the coupling in Its trace 
is visible in the explicit structure of all the diagonal two- loop anomalous dimensions. More- 
over, this argument naturally emerges when using dispersive reasoning to carry out a careful 
treatment of the appearance of the running coupling in inclusive processes [10]. 

We are still far from a good understanding of how to simplify the structure of multi- 
loop anomalous dimensions, notably because of complications that arise from off-diagonal 
transitions (see appendix). Nevertheless the belief that the x — >■ 1 limit is under control 
has led us to investigate the implications of eq. © for the large- structure of three loop 
anomalous dimensions. 

As we shall see this will provide insight concerning a "very suggestive" relation noted by 
Moch, Vermaseren and Vogt (MVV) which, in their words, ^^seems to call for a structural 
explanation": in both the non-singlet quark [2] and diagonal singlet quark-quark and gluon- 
gluon three-loop splitting functions [3], they observed that the third-loop coefficients, Cf, in 

^Such a reformulation of the notion of parton splitting functions originally proposed in [7] has been carried out in 
detail in the context of heavy-quark fragmentation functions, where it was found to greatly improve the perturbative 
series [8]. 




(3) 



where a = —1 (-1-1) for the space-like (time-like) case. The choice of z'^Q'^ as the logarithmic 
ordering parameter (parton evolution 'time') corresponds to ordering parton splittings in the 
fluctuation lifetimes of successive virtual parton states. Assuming the new 'P(x) splitting 
functions to be identical in the space-like and time-like cases one obtains that the 'traditional' 
two-loop splitting functions in the two cases, 

p(2)-r(x) differ by 
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the large- A'^ expansion of the n-loop anomalous dimensions,^ 



7aa{N) = -A"" {In N + -fe) + - N~HnN + 0{N-^) , 

7..-E7n.»(^)"' ^-E<(^)^ etc. (5) 

are simply related with A2 and Af, i.e. Cf = 2AfA2, where a = q,g. This supplements the 
two-loop relation, Cf = (^5^)^ [4]. 



2 MVV relation 

For the purpose of studying the x — > 1 limit we initially approximate 'P{x, as) by the product of 
the physical coupling, Oph, and the 1-loop splitting function, V{x,as) ^ (aph/47r)P(^) (x) (for 
compactness we will write (aph/47r) = a). To deal with the correlated z and dependences 
in the right-hand side of ©, we rewrite it as 

dtD{x,Q^) = [' —P^^\z) [e''^-^(-)^-al {e'^'^'^W (-,Q^)\ , (6) 
Jo z L J L \z J i 

where both a and D are now evaluated at scale and /?(a) = —da/dt = — l^„=o /^"'^""''^ • 
The Mellin transform of this equation results in the formal expression, 

dtDM = 7iiN + Pia)d^ + adt) aD^ , (7) 

where 7i(A^) is the Mellin transform of the first order splitting function P^^\x). We note 
that dt operates only on and not on a. This gives an all-order model for the anomalous 
dimension, 7(A^) = Dj^^dtD]^. 
Expanding ((7| results in 

7EE7[a] = a-^,+^^D-^{l3da, + adt){aD) + +\^^D-\pda + adtf{aD) + ... 
= 071 + 7i(/3 + c^a7) + \li D'^iPda + crdt){/3D + aaDj) + ... 
= a7i+7i(/3 + aQ7) + i7i [072 + ^(2^7 + Q/35„7) + /?9«/3] +0(a^), (8) 

where dots indicate derivatives with respect to and 71 = 71 (A^). Solving this iteratively 
produces 

7 = a7i+a2 j^{f3o+aj,)+a^ [j,{f3i + a^,{l3o + aj,)) + ^71(7^ + Sa/Soji + 2/3^)] +0{a^) . 

(9) 

For the purpose of understanding the MVV relation it suffices to take 71 = —Ai In + 0{1) 
and to keep in © only the term oc 7171, giving 

7 = -aAi \nN + const. + aa^Aj + 0{N-^) . (10) 



^We define jn{N) = dzz^ '^pM (2); this has the opposite sign to the convention of MVV; we have also changed 
the sign in front as compared to eq. (3.10) of ref. [2], which contains a misprint. 
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Recalling that A^a = A'' (with Al = ACp and Af = 4Ca), we can then write the followinK 
all-order relation between C"^ and A^, 

C^ = -a{A^f, (11a) 

or equivalently, in terms of the expansion coefficients Cn, 

Ci = 0, C2 = -aAj , Cs = -2aAiA2 , C4 = -t{AI + 2AiA^) , etc. (lib) 

where we have suppressed the index a = q,g. For the space-like case {a = —1) this explains 
the MVV observation. For the time-like case we have only the two-loop result [4] to compare 
to, and it agrees. 



3 Pushing our luck 

Motivated by the idea that the universality of soft gluon emission holds both in singular and 
constant terms in gluon energy [9], one may attempt to trace further terms of the large- A'" 
expansion generated by eq. ©• We definitely expect this push to fail at the level of 
(possibly modulo logarithms, see below) because this corresponds to 'quantum' terms in the 
splitting function, which vanish as 1 — x. However we would expect to have control over the 
1/N term in the anomalous dimension, 

7(iV) = -Aii^iN+l) + 7e) + -B - CiipiN) + 7e) N'^ + D N'^ + 0{N-^ log^ N) . (12) 

Compared to © we have shifted the argument of the logarithm in the A-term, N + 1, 
added the constant 7e in the C-term and then replaced logarithms with ip functions. These 
modifications do not affect the first three functions A, B and C but serve to simplify the next 
subleading term oc 1/A^. Additionally they lead to a compact x-space image of (jl2p . 

p(x) = , +B5{l-x) +Cln(l-x) + D + 0((l-x)logP(l -x)) . (13) 

(l-x)+ 

Here, the presence of x/(l — x) in the first term (as opposed to 1/(1 — x)) is a consequence of 
Low's theorem [9]. In the calculation of the C coefficients we could safely ignore the 0(1) piece 
of 71. This is no longer possible when calculating D because, in our non-linear construction, 
this constant {unity in Mellin space) is multiplied by 71 ~ 1/A^ thus contributing to D. The 
extension of Q to account for B in all orders is obtained by generalising 

OTi ^ a(7i -Bi) + B = -^(ln(A^+l) + 7e) + 5 + 0{N~^) . (14) 

The reason why the structure of the Taylor expansion Q is unmodified modulo this simple 
substitution is that B, being a constant, disappears everywhere but un-dotted factors of 7j. 
This leads to the following all-order expectation for D, 

= - <yB-^ , (15) 

where we have rewritten the A'fP term that comes from eq. (jH)) as —{dtA'^) (recalling the 
definition of /? in terms of the physical coupling). The MS expansion for D"" is then 

1)1 = 0, D2 = -AiiaBi + po) , D3 = -Ai{aB2+Pi)-A2iaBi + 2-po), etc. (16) 
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Hard luck. Examining the full known results for the two and three-loop splitting functions, 
we find agreement for D2 (space and time-like, and quark and gluon channels); however the 
result for the space-like -D3 is as follows (for both quarks and gluons) 

D^ = Ai{B2-Pi) + A2{Bi-l-(3o). (17) 

There is one mismatch between eqs. ()16p and (|17p . which we have highlighted in boldface. Had 
((T3|) contained (9faMs)/aMs instead of {dtA°')/A'^ we would have obtained agreement with the 
full result for D3, however we see no reason why it should be the MS coupling that appears 
there instead of the physical coupling (which is equivalent to putting A""). The remarkable 
simplicity of the mismatch calls for a further structural explanation. 

Good luck. Despite the disagreement in the comparison with the exactly calculated sub- 
leading D term (which we hope can be understood) we have also investigated the coefficient 
of terms that vanish for a; — > 1 but that are logarithmically enhanced there. We have found 
agreement using Q for the oi1{\ — x)ln(l — x) and the 0^(1 — x)ln^(l — x) terms in the 
(space-like) non-singlet anomalous dimensions, while a^(l — x)ln(l — x) contains structures 
that remain to be understood. As for the diagonal singlet anomalous dimensions, at three loops 
the coefficient of a\{l — x) ln^(l — x) agrees only in its n ^-independent parts and there is an 
additional unexpected a^(l — x) ln^(l — x) contribution proportional to nj, whose origin may 
only be explained once the higher-order structure of the off-diagonal splittings is elucidated. 

The wealth of structure that is present in higher-order splitting functions is suggestive of 
underlying simplicity. Possible sources of such simplicity, as proposed here, are the universal 
nature of soft gluon radiation and the reformulation of the notion of parton splitting functions 
with the aim of preserving universality between space and time-like parton multiplication 
(Gribov-Lipatov reciprocity). Whether this picture can be made fully consistent remains to 
be seen. We look forward to future work shedding more light on this question. 
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Appendix 

We do not know how to generalise (jSJ to non-diagonal transitions. These are considerably more 
divergent at large x [13] than would be expected based on the non-linear relations that we 
propose here, going as a" ln^"~^(l— x) at large x. We suspect that the origin of these additional 
logarithms may be that in the MS factorisation scheme, for a — > 6 transitions with a ^ b, 
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the splitting functions could pick up residues from ratios of non-cancelling divergent Sudakov 
exponents (as well as from singular integrals of these ratios). This belief is not inconsistent with 
the MVV observation that in the supersymmetric case most of these logarithmic enhancements 
cancel, since in this case the Sudakov exponents become identical for quarks and gluons and 
do cancel. 

These problems of non-diagonal terms may be responsible for the following fact: at two 
loops, the V that appears in (j^J for gluon-gluon splitting is universal (identical for space and 
time-like cases) only for two of the colour structures, C\ and CAnf. The analogue of Q for 
the remaining colour structure, CpUf relates gluon-gluon and singlet quark-quark splittings. 

On the left-hand side one finds the combinations Pgg^'^ — Pqq^'^ and Pqq^''^ — P^'^ , while 
on the right-hand side one has convolutions involving Pg]j^{x/z) ■ InzPgg^(z) and Pj^j^{x/z) ■ 
\nzPj,^{z) [7]. 
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